Metal-Insulator transition in Hubbard-like models with random hopping 
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An instability of a diffusive Fermi liquid, indicative of a metal-insulator transition (expected 
to be of first order), arising solely from the competition between quenched disorder and short- 
ranged interparticle interactions is identified in Hubbard-like models for spinless fermions, subject 
to (complex) random hopping at half-filling on bipartite lattices. The instability, found within a 
Finkel'stein non-linear sigma model treatment in d — (2 + e) > 2 dimensions, originates from an 
underlying particle-hole like (so-called chiral) symmetry, shared by both disorder and interactions. 
In the clean, interacting Fermi liquid this symmetry is responsible for the (completely different) 
nesting instability. 

PACS numbers: 71.30.+h, 71.10.Fd, 72.15.Rn 
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Understanding the combined effects of quenched dis- 
order and interparticle interactions in electronic systems 
remains one of the central problems in solid state physics. 
Models of noninteracting electrons subject to static, ran- 
dom impurity potentials provide the simplest descrip- 
tion of disordered metals; many analytical and numerical 
studies have shown that such models exhibit a continu- 
ous metal-insulator transition (MIT) in three dimensions 
(3D)ji accessed by varying either the disorder strength, 
or the Fermi energy relative to the mobility edge. By 
contrast, all electronic states are typically exponentially 
localized in one and two dimensions by arbitrarily weak 
disorder. (Notable exceptions occur in the presence of 
spin-orbit scattering, and in the systems discussed, e.g., 
in Refs.ii.) 

The above description ignores the important effects of 
interparticle interactions. Unfortunately, however, theo- 
ries capturing the competition of both disorder and inter- 
actions are typically quite complex, and difficult to ana- 
lyze reliably^ These issues have come again to the fore- 
front of scientific debate in view of discussions concerning 
a metal-insulator transition in 2D, and related experi- 
mental results on 2D semiconductor inversion layers. 6 

In the present work, we identify a novel "Anderson- 
Mott" instability of a diffusive Fermi liquid in 
d = (2 + e) > 2 spatial dimensions (e -C 1), which 
arises solely from the competition between disorder and 
short-range interactions. This instability is indicative of 
a metal-insulator transition (expected to be of first or- 
der) from the diffusive Fermi liquid to an insulating state 
dominated by both strong disorder and interactions. (See 
Fig- HI below.) Since the system that we study has no lo- 
calized phase with disorder in the absence of interactions, 
a localized phase can only appear due to the presence of 
the interactions. We expect our result to be relevant for 
sufficiently strong disorder in d — 3 dimensions. 

Specifically, we analyze a class of "Hubbard-like" 
models^ of spinless fermions, at half- filling on bipartite 
lattices, with random (short-range) hopping between the 
two sublattices. [The Hamiltonian is given below in 
Eq. ([I}.] In every realization of the (static) disorder, such 
a model possesses a special particle-hole like symmetry, 



which we will refer to as sublattice symmetry (SLS). [SLS 
is termed "chiral symmetry" in the classification scheme 
of Ref. | (see also Refs. |^laf9|JTo|JTllIT^fT^) .] 

In the absence of disorder, it is SLS which is respon- 
sible for the "nesting" condition of the Fermi surface. 
Fermi surface nesting is, in a sense, the defining prop- 
erty of (clean) Hubbard-like models for interacting lattice 
fermions in d > 2. It is the nesting condition which makes 
the ballistic Fermi liquid phase at half filling in such mod- 
els unstable to Mott insulating order in the presence of 
generic, arbitrarily weak interparticle interactions ^ 4 ' 15 : 16 
Here we consider, in addition, complex random (nearest- 
neighbor) hopping, which breaks time reversal invariance 
(TRI) in every realization of disorder^ For our system 
of spinless fermions, this is consistent with the applica- 
tion of a random magnetic field to the otherwise clean 
model. 

Our principal motivation for studying such a model is 
that we expect both disorder and interparticle interac- 
tions to play important roles in the description of the 
low-energy physics. Because random hopping preserves 
the special SLS, our disordered model retains the nesting 
instability of the associated clean system. This instabil- 
ity can therefore compete with the unusual localization 
physics of the disordered, but noninteracting model (sec 
below) . The further assumption of broken TRI guaran- 
tees that we do not have to confront an additional super- 
conducting instability ^ 14 ' 15 ! 16 We expect the instability 
that we identify in this work in the simultaneous pres- 
ence of both disorder and interactions to occur in three 
dimensions for sufficiently strong disorder, and we stress 
that it is clearly distinct from the pure Mott nesting in- 
stability, although the latter also appears in our model 
phase diagram (see Fig. [1] below). We note that the ef- 
fects of hopping disorder upon the Neel ground state of 
the (slightly more complex) spin-1/2 Hubbard model at 
half filling were studied numerically in Ref. Ila . although 
these studies were limited to d = 2. 

A second motivating factor is that, interestingly, the 
presence of SLS radically changes the localization physics 
of the disordered, noninteracting random hopping model 
[Eq. {TJ, below, with V = U = 0]. SLS enables the ran- 
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dom hopping model to evade the phenomenon of Ander- 
son localization. Specifically, the noninteracting system 
exhibits a critical, delocalized phase at the band center 
(half filling) in one, two, and three dimensions for finite 
disorder strength, with a strongly divergent low-energy 
density of states in d — 1, 2.^ & 10 ' 11 ' 12 In particular, 
there is no MIT and no Anderson insulating phase in 
d = 3 (in the absence of interactions). Random hopping 
models have been of significant theoretical interest in the 
recent past, both because of the unusual derealization 
physics described above, but also because these models 
have proven amenable to a variety of powerful analyti- 
cal techniques in d < 2, with many exact and/or non- 
perturbative features now understood ^ 10 ? 12 This situ- 
ation should be contrasted with our understanding of 
the conventional noninteracting ( "Wigner-Dyson" ) MIT, 
which is based largely on perturbative results in d > 2. 1 
Our work here addresses the effects of interparticle inter- 
actions in random hopping models for d > 2. 

Our starting point is the following extended Hubbard- 
like Hamiltonian for spinless fermions hopping on a bi- 
partite lattice at half filling: 



Ai C Bj 



H.c. + 2V^y^ SfiAiSfiBj 

(ij) (ij) 

+ U C^2 ^AiSfiAi' + ^ $nBj&n B j>)- (1) 

<(»'» (W)) 

Any bipartite lattice may be divided into two interpene- 
trating sublattices, which we distinguish with the labels 
A and B. In Eq. ([T|), c' Ai and c B - denote fermion creation 
and annihilation operators on the A and B sublattices, 
respectively. Here, i and j respectively index the A and B 
sublattice sites, and the sums on (ij) run over all nearest 
neighbor A-B lattice bonds, while the sums on ((ii 1 )) and 
((jj')) run over all next-nearest neighbor (same sublat- 
tice) pairs of sites. The homogeneous hopping amplitude 
t in Eq. ([T]) is taken to be purely real; disorder appears 
in the perturbation Stij. We take the amplitude Stij to 
be a Gaussian complex random variable with zero mean, 
statistically independent on different lattice links. The 
operators Stia/b = ( c a/b c a/b ~~ \) m © denote the 
deviations of the local sublattice fermion densities from 
their value at half filling; the interaction strengths V and 
U appearing in this equation couple to nearest neighbor 
and next-nearest neighbor density-density interactions, 
respectively. 

The Hamiltonian H in Eq. |T]) is invariant under the 
(antiunitary) sublattice symmetry (SLS) transformation 



^Ai ' C B j 



(2) 



(all complex scalar terms in H are complex 
conjugated)*^ In the clean limit, Eq. |T|) 
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for all lattice bonds (ij), SLS is responsible for the 
nesting condition of the noninteracting Fermi surface. 
As a result of nesting, the Fermi liquid phase of the 
clean model is unstable to charge density wave (CDW) 
order for any 2V > U > 0^2. 



The effects of the interparticle interactions U and 
V upon the delocalized phase of the disordered 
model given by Eq. |T]) may be investigated by us- 
ing Finkel'stein's generalized non-linear sigma model 
(FNLctM) approach^ to formulate the low-energy effec- 
tive continuum field theory. The latter can be studied 
using a controlled e expansion in d = 2 + e dimensions 
(0 < e < 1). We use the Schwinger-Keldysh 21 method 
to perform the ensemble average over realizations of the 
hopping disorder. The FNLctM is derived following the 
standard methodology^^ 

In the present case, the resulting FNLctM is described 
by the generating functional 2 ^ 



Z = \VQe 



-Si 



(3) 



where 
S D = 



/ d d rTr 
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VQt( r ). VQ(r) 
+h j rf d rTr{^[Qt(r) + Q(r)]} 
-g^j / rf d r{Tr[Qt (r)V g(r)]} 



(4) 



and 



^ = l E Z a I dtd d r{2T v Qtr(r)Qtt(r) 

a=l,2 ^ 

+r £/ [Q-(r)Q-(r) + Qjf 1 WlrW]} 



The field variable 



Q(r)-Q#(r) 



(5) 



(6) 



in Eqs. ©-(O is a complex, "infinite-dimensional" 
square matrix living in d spatial dimensions, where in- 
dices t and t' belong to a continuous time or (via Fourier 
transform) frequency space, and where indices a and b 
belong to a 2-dimensional "Keldysh" species space, with 
a,b E {1,2}. 21 ' 22 In Eq. (g]), Tr denotes a matrix trace 
over time (or frequency) and Keldysh indices. Q(r) sat- 
isfies in addition the unitary constraint 



Qt(r)Q(r) = 1. 



(7) 



The matrix Q and its adjoint may be interpreted 2 ^ as 
continuum versions of the same- sublattice fermion bilin- 



ears 



Q1^c A (t)c^(t'), 
Q\v - c%(t)c%(t'). 



(8a) 
(8b) 



The action given by Eq. ([J| describes the low-energy 
diffusive physics of the noninteracting random hopping 
model; a replica version of the noninteracting sigma 
model with action was originally studied by Gade 
and Wegner. 3 This noninteracting sector of the FNLctM 
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contains three coupling constants: A, Xa, and h. The pa- 
rameter 1/A is proportional to the dimcnsionlcss dc con- 
ductance g of the system (i.e., is related to the disorder 
strength), while Xa denotes a second measure of disorder, 
unique to this sublattice symmetry class, which strongly 
influences the single-particle density of states^ ' 10 ' 11 The 
parameter A^i may be simply interpreted as character- 
izing the strength of long-wavelength, quenched random 
orientational fluctuations in bond strength dimerization 
of the random hopping model. 20 Finally, h is a dynamic 
scale factor, which determines the dynamical critical ex- 
ponent z in Eq. ()10c|) . below, through the condition 
dlnh/dl = 0i£ i 13 i 20 The interparticle interactions appear 
in Si, defined by Eq. ((5|). Given Eq. (J3]), we may interpret 
Qtt(r) and Q} t aa (r) as continuum local density operators 
on the A and B sublattices, respectively. Then the inter- 
action couplings IV and Tjj in Eq. ((5]) describe generic 
short-ranged intersublattice and same-sublattice density- 
density interactions, respectively [compare to Eq. |T]), 
above]. Finally, £ a = ±1 in Eq. ([5]), for Keldysh species 
label a = 1, 2. 

Using a Wilsonian frequency-momentum shell back- 
ground field methodology^ we have performed a one- 
loop renormalization group calculation on the model 
defined by Eqs. (J3])-©. The calculation is straight- 
forward, though rather lengthy; the details will be pub- 
lished elsewhere. 20 Below we simply state our results. In 
order to do so, it is convenient to introduce the following 
effective interaction couplings 



7s 



v), 7c 



(9) 



The interaction strength 7 S couples to the square of the 
(smooth) local charge density in the continuum theory, 
while 7 C couples to the square of the sublattice staggered 
charge density. In accordance with the discussion in the 
paragraph below Eq. ([21), we expect j c < to promote 
charge density wave (CDW) formation, while j c > 
should suppress it. 

We find the following one-loop RG flow equations for 
the couplings A, A^, 7 S , 7 C , and h in d = (2 + e) dimen- 
sions: 



dX 

m 



IT = -eA-A 2 7c 



-2A 2 



1 



7s 



ln(l- 7 s) 



dX A 
dl 

dfh 
di 

drfc 
dl 

d\nh 
dl 



, 9 „ A a dX 
= eA^ + A +2— — , 

= Aa(1 - 7s) (7s + 2 7c - 27 s 7 C ) 
-A(l - 7s ) ( 7s + 7c - 27 s 7c) , 

= X A (7c + 2 7s ) - A ( 7s + 7c) 
+A[2 7 cln(l- 7s )+ 7s 7c] - 2 7 2 , 

= (d- z) + X A + X(-f c -j s ). 



(10a) 
(10b) 

(10c) 

(lOd) 
(lOe) 



Here, I is the logarithm of the spatial length scale. These 
flow equations are given to the lowest non-trivial order 
in the couplings A, A^, and 7 C , but contain contributions 
from 7 S to all orders; FinkcPstcin's NLerM formulation 
provides* a perturbative expansion which is controlled 
by the (small) dimensionless resistance A, but does not 
require the interaction strength 7 S to be small. Before 
turning to an analysis of our results, Eqs. (|10ap - (|10ep . 
we provide interpretations for various key terms appear- 
ing in them. First, the term in square brackets on the 
second line of Eq. (|10a|) is the usual correction to the di- 
mensionless dc resistance A, arising from the short-ranged 
interparticle interactions ^ 23 ' 24 and corresponds^ to co- 
herent backscattering of carriers off of disorder-induced 
Friedel oscillations in the background electronic charge 
density^ The last term in Eq. (|10d|) drives the CDW 
instability, which is a remnant of the clean Hubbard-like 
model (recall that in our conventions 7 C < signals this 
instability) . 

Now we analyze our results. In d = 2 dimensions, 
integrating Eq. (|10p for generic initial conditions shows 
that the critical, delocalized phase of the half-filled, non- 
interacting random hopping model 3 - is unstable to the 
effects of short-ranged interparticle interactions. We find 
that either 7 C — > — oo signaling CDW formation, or that 
A, Xa — * oo and 7 C — » +oo, indicating a flow toward both 
strong disorder and strong interactions. Regardless, we 
expect the 2D interacting, disordered Hubbard model to 
be an insulator at zero temperature. This should be com- 
pared to an analogous result 1 - 3 - previously obtained for a 
TRI, interacting random hopping model on the honey- 
comb lattice. This physics is consistent with numeri- 
cal studies 1 ^ of the half-filled spin- 1/2 Hubbard model in 
d = 2, which have shown that TRI random hopping dis- 
order preserves the compressibility gap of the clean Mott 
insulator, and that the disordered and interacting system 
shows no signs of metallic behavior. 

The situation in d = (2 + e) > 2 dimensions is more in- 
teresting. Upon increasing e from zero, a narrow, irregu- 
larly shaped sliver corresponding to a stable metallic, dif- 
fusive Fermi liquid state opens up in the four-dimensional 
(A, Aa, 7s, 7c) coupling constant space. The sliver en- 
closes the line A = Xa — 7c = 0, with — oo < 7 S < 1, the 
entirety of which is perturbatively accessible because the 
FNLctM does not require the interaction strength 7 S to 
be small. A highly schematic 3D "projected" phase dia- 
gram is depicted in Fig. [TJ In this figure, the interaction 
constants reside in the horizontal plane, while the ver- 
tical direction schematically represents (both) disorder 
strengths; the shaded sheath is a cartoon for the bound- 
ary of the stable metallic region, which resides between 
it and the ballistic (A = Xa — 0) plane. The "height" 
of the stable metallic region in the "disorder" directions 
(A, Xa) is controlled by e, although the precise shape and 
size of the phase boundary varies with 7 S , and is difficult 
to characterize analytically. Over the range of perturba- 
tively small values of j c , the stable Fermi liquid phase 
resides in the region 7 C > 0, and terminates near 7 C = 0. 
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FIG. 1: Schematic phase diagram in d = 2+e, with < e < 1. 
A and Aa are measures of the disorder, whereas 7„ and 7 C 
characterize the interaction strengths. The stable metallic 
phase resides between the ballistic plane (A = \a — 0) and the 
shaded sheath; the thick arrows indicate the two instabilities 
of the (metallic) diffusive Fermi liquid discussed in the text. 

The flow equations ()10a|) - (|10dj) possess no perturba- 
tively accessible, nontrivial RG fixed points for d > 2, 
and thus no continuous metal-insulator transition can be 
identified. However, the two instabilities described above 
for the 2D case persist for d > 2, and become clearly dis- 
tinct roads out of the metallic state. The conventional 
CDW instability always occurs for initial j c < and suf- 
ficiently weak disorder, i.e., when A, A^ <C e, and is rep- 
resented by the flow j c — > — oo. This flow is accompanied 
by a decay in both disorder strengths A, A^. 



The primary result of this paper is the identification 
of a second route out of the diffusive Fermi liquid phase 
in d — (2 + e) > 2 dimensions, independent of the 
Mott CDW instability, arising solely from the competi- 
tion of disorder and interaction effects. As in the 2D 
case, this second route is characterized by a flow off 
to both strong disorder (A, A^ — > oo) and strong inter- 
actions (7 C — ► +oo), as indicated by the thick arrows 
emerging from the 7 C > portion of the phase boundary 
shown in Fig. [TJ we call it an Anderson-Mott instability. 
Even though there is no perturbatively accessible fixed 
point, this new instability is nonetheless perturbatively 
controlled in d = (2 + e) over a wide range of initial con- 
ditions when e -C 1; in particular, it is perturbatively ac- 
cessible over the entire range < % < li^Z Numerically 
integrating Eqs. (|10a|) - (|10d[) for small e < 1, we find that 
the Anderson-Mott instability can apparently always be 
reached by increasing only the dimensionless resistance A. 
We expect the boundary separating the flow toward the 
stable metallic regime from that toward the regime of the 
Anderson-Mott instability to represent a disorder-driven, 
first order metal- insulator transition (MIT). We empha- 
size that a MIT does not exist in the noninteracting ran- 
dom hopping model, which possesses only a delocalized 
phase at half-filling for finite disorder in d > 1,— £ i 10 i 1:L 
while the clean spinless Hubbard model possesses only 
the Mott CDW instability. RG flow equations in related 
systems of spin-1/2 fermions were recently obtained in- 
dependently in Ref. \2Q. 
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